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PACS 05 . 10 . -a - Computational methods in statistical physics and nonlinear dynamics 

PACS 87. 15. A — Theory, modeling, and computer simulation 

PACS 87. 15. Cc - Folding: thermodynamics, statistical mechanics, models, and pathways 

Abstract. - By means of multicanonical computer simulations, we investigate thermodynamic 
properties of the aggregation of interacting semiflexible polymers. We analyze a mesoscopic bead- 
stick model, where nonbonded monomers interact via Lennard- Jones forces. Aggregation turns 
out to be a process, in which the constituents experience strong structural fluctuations, similar to 
peptides in coupled folding-binding cluster formation processes. In contrast to a recently studied 
related proteinlike hydrophobic-polar heteropolymer model, aggregation and crystallization are 
separate processes for a homopolymer with the same small bending rigidity. Rather stiff semiflex- 
ible polymers form a liquid- crystal- like phase, as expected. In analogy to the heteropolymer study, 
we find that the first-order-like aggregation transition of the complexes is accompanied by strong 
system-size dependent hierarchical surface effects. In consequence, the polymer aggregation is a 
phase-separation process with entropy reduction. 
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Cluster formation and crystallization of polymers are 
processes which are interesting for technological applica- 
tions, e.g., for the design of new materials with certain me- 
chanical properties or nanoelectronic organic devices and 
polymeric solar cells. From a biophysical point of view, 
the understanding of peptide oligomerization, but also the 
(de)fragmentation in semiflexible biopolymer systems like 
actin networks is of substantial relevance. This requires a 
systematic analysis of the basic properties of the polymeric 
cluster formation processes, in particular, for small poly- 
mer complexes on the nanoscalc, where surface effects arc 
competing noticeably with structure-formation processes 
in the interior of the aggregate. 

A further nrotivation for investigating the aggrega- 
tion transition of semiflexible Ziomopolymer chains derives 
from the intriguing results of our recent study of a similar 
aggregation process for peptides [TJH], which were mod- 
eled as /seieropolymers with a sequence of two types of 
monomers, hydrophobic (A) and hydrophilic ones (B). By 



M E-mail 
( b ) E-mail 
M E-mail 



junghans@mpip-mainz . mpg . de 
m.bachmann@fz-juelich.de 
Wolfhard . JankeSitp . uni-leipzig . de 



specializing the previously employed heteropolymer model 
to the apparently simpler homopolymer case, we aim by 
comparison at isolating those properties which are mainly 
driven by the sequence properties of heteropolymers. In 
fact, while in both cases the aggregation transition is a 
phase-separation process, we will show below that for ho- 
mopolymers the aggregation and crystallization (if any) 
are separate conformational transitions - unlike our study 
of heteropolymer aggregates where they were found to 
coincide [UH]. The physical origin causing these differ- 
ences will be explained within the microcanonical formal- 
ism 3,4 , which proves [IJ[2] to be particularly suitable for 
this type of problem. 

We thus consider the same model as in PQ[5], but here 
we assume that all monomers i^ = 1, . . . , NU*) of the /zth 
chain (fi = 1,...,M) at positions x,- are hydrophobic 
(A). The bonds between adjacent monomers are taken 
to be rigid (bead-stick model) and pairwise interactions 
among nonbonded monomers are modeled by a Lennard- 
Jones potential 
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where r^ j u = |x* —Xj J is the distance between monomers 
i^ and j v of the /xth and i/th chain, respectively. In the 
peptide model, the coefficient of the van der Waals con- 
tribution oc r~ 6 would depend on the type of monomers 
involved [1,2 . Intra-chain (/x = u) and inter-chain (/i ^ v) 
contacts are not distinguished energetically. The scmiflex- 
ibility of a chain is described by the bending energy 



41d = «E( 1 - cos ^)' 
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where < $j l < ir is the bending angle formed by the 
monomers i M , i^ + 1, and i^ + 2. For the comparison 
with our recent heteropolymer aggregation studies OH], 
we choose in most simulations a bending rigidity k = 0.25, 
which is at the rather floppy end of semiflcxibihty. Thus, 
the single-chain energy reads 



E M = R 
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and the total energy of the polymer system is given by 



b-E bM + EE^(»u) 



(4) 
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All chains are assumed to have the same degree of poly- 
merisation, i.e., the same number of monomers, JVM — TV, 
H = 1,...,M. 

We have performed multicanonical computer simula- 
tions [5] for different system sizes. The multicanonical 
weights were calculated recursively [5] in 50 iterations with 
5 x I0 6 x M 2 single updates (spherical pivot rotations [7], 
semilocal crankshaft moves ) each. The production run 
with fixed multicanonical weights included 5 x I 8 x M 2 
single updates. The generic result obtained from these 
simulations is a precise estimate for the density of states 
g(E) or microcanonical entropy S(E) = ks \ng(E), which 
will be central in the second part of our data analysis. 

From a more standard canonical perspective, confor- 
mational transitions between structural macrostates are 
signalized by peaks in the temperature-dependent fluctu- 
ations of energy, i.e., the specific heat per monomer 



c v = {(E 2 )-{E) 2 )/N tot k B T, 
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where iVtot = NM and ks = I in our units. Know- 
ing g(E), this can be straightforwardly computed for any 
temperature T. In Fig. 02a), the specific-heat curve for a 
system of two identical semiflcxiblc polymers (2 x ^13) 
is compared with the energetic fluctuations of a single 
chain (I x ^13). The single chain exhibits a very weak 
coil-globule collapse transition (shoulder near T rs 0.88), 
whereas the crystallization near T w 0.24 is a pronounced, 
separate process. The thermodynamic phase behavior of 
single semiflexible polymers in solvent has already been 
subject of numerous studies, with particular focus on stiff- 
ness and finite chain length effects [5HTU]. where it was 
shown that the globule-solid transition is more influenced 



by stiffness effects than the coil-globule transition. In 
particular, for longer chains, it was found that, depend- 
ing on the stiffness, single collapsed semiflexible polymers 
form spherical, ellipsoidal, and disklike globules, as well 
as toroids |10) . Our first result for the semiflexible multi- 
ple-chain system read off from FigOJa) is that aggregation 
and collapse are not separate processes (near T w 0.97), 
similar to the corresponding heteropolymer system 00 . 

At about T en 0.24 (close to the single-chain freezing 
temperature), the multiple-chain homopolymer complex 
crystallizes in a separate process. This is in strong con- 
trast to the heteropolymer systems, where aggregation, 
collapse, and crystallization (hydrophobic-core formation) 
is a single-step process pQ[2] • In Fig. O a ) , we have also 
included a comparison with a two-chain system with much 
stronger bending rigidity k = 10. As expected, there is a 
single aggregation transition near the temperature, where 
the system of less stiff semiflexible polymers with k — 0.25 
collapses. However, a further crystallization process at 
lower temperatures does not occur: There is no globular 
(liquid) pseudophase of defragmented relatively stiff semi- 
flexible polymers. 

Of course, in finite and small polymer systems confor- 
mational macrostates form no "phases" in a strict ther- 
modynamic sense. We hence call the stable macrostates 
in the following "pseudophases" in order to emphasize the 
difference. Nonetheless, by introducing a "phase" separa- 
tion parameter 
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where r^Q M = ^2" , r^ /N [2] is the center of mass of 
the polymers, a clear discrimination of the different pseu- 
dophases can be made. Small values of T correspond to 
aggregated and higher values to fragmented conforma- 
tions. In Fig. (Ub), the canonical expectation value (r) 
and its fluctuation d(T)/dT are plotted. The peak posi- 
tion of d{T)/dT coincides nicely with the corresponding 
peak temperature of the specific heat and thus signals the 
aggregation transition. 

In Fig. [21 the aggregation parameter and its fluc- 
tuations are shown for different system sizes of up to 
four chains. The peak shifts towards higher tempera- 
tures and gets sharper with increasing number of chains, 
since intra-chain and inter-chain monomer-monomer con- 
tacts are not energetically distinguished. The hypothetic 
maximal total number of intrinsic contacts is n™t? a = 
M(N - 2)(N - l)/2 - MN 2 = NN tot , and for the max- 
imally possible number of inter-chain contacts one has 
n m j* c x r = M(M - l)N 2 /2 ~ M 2 N 2 = Af t 2 ot [TTJ. For large 
M, the relative fraction renter of the inter-chain contacts, 
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2xAi 3 , k = 10 

2xA 13 , k = 0.25 
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Fig. 1: (a) Specific heat per monomer as a function of temper- 
ature for a single semiflexible homopolymer with 13 monomers 
(1 x Ais, k = 0.25) and for systems of two such chains (2 x A13) 
with different bending rigidities, (b) Canonical expectation 
value (r) and fluctuation d(T)/dT of the aggregation parameter 
F, defined in Eq. (|6j), for the two-chain system with k = 0.25. 



behaves like ri nte r ~ 1 — A/ -1 , i.e., the relative influence of 
inter-chain contacts increases rapidly towards unity with 
the number of chains. In consequence, aggregation domi- 
nates over collapse of the individual chains. Even for the 
two-chain system 2 x A13 this estimate is reasonable: The 
energy of the lowest-energy conformation we found nu- 
merically is £'( mm ) « —83.61 and the contribution of the 
inter-chain contacts is E, 



with 



(min) 



( min ) ~ cnon - ~( min ) p(min) 
inter bU - ZU ~ 7 inter ^ 

0.60. This coincides nicely with the corre- 
sponding value of the above contact ratio, ri n t er ~ 0.56. 
In fact, considering energetic and structural fluctuations 
of the larger systems with three and four chains, we have 
not found indications for an additional collapse transition 
at temperatures higher than the aggregation transition. 

Below the aggregation temperature, the entropic loss 
of the individual chains is overcompensated by the ener- 
getic gain of forming a joint globular aggregate. However, 
the entropic change while passing the aggregation transi- 
tion is noticeably smaller than what we found recently for 
heteropolymer systems, where no separate freezing transi- 
tion occurs P]. In the intermediate fluid "globular" pseu- 
dophase, the aggregate of the homopolymers thus behaves 
like a single chain of length MN. Consequently, reducing 
the tenrpcrature, the aggregate optimizes the monomer ar- 
rangements in order to maximize energetic contacts. In- 
dicated by the peak in the specific heat near T m 0.2 
(see Fig. [T] and also Ref. [7]), the small globular aggre- 
gates freeze into spherical amorphous structures with a 
maximum number of inter-chain contacts. For rather stiff 
semiflexible polymers (as our example with k = 10), how- 
ever, a separate freezing transition does not occur and the 
peak of cy near T ps 0.9 indicates a single-step transition 
from rod- like coils to a liquid-crystal-like phase |10j . 




Fig. 2: (a) Canonical expectation values and (b) fluctuations 
of the aggregation parameter F for multi-chain systems with 
K = 0.25. In (a), also an exemplified globular 4 x A13 aggregate 
is depicted. 



In our recent studies of heteropolymer aggregation [H2] , 
we also observed a pronounced single transition of a differ- 
ent nature. In this case, the formation of a heteropolymer 
complex consisting of different chains with compact hy- 
drophobic core is roughly a single-step process, because 
hydrophobic-core formation ( "freezing" ) favors conforma- 
tions with very small entropy. For semiflexible homopoly- 
mer systems, we find in the rather floppy limit that the 
freezing temperature (T w 0.2, for single chains or globu- 
lar aggregates with k = 0.25) is almost identical with the 
aggregation tenrpcrature for heteropolymer systems [U[2]. 
This coincidence in the behavior of these different sys- 
tems is due to the formation of a single, very compact 
hydrophobic domain (the monomers of the interacting ho- 
mopolymers are as hydophobic as the A-type monomers 
of the heteropolymer) which maximizes the number of en- 
ergetic contacts. Stiff homopolymers, on the other hand, 
cannot form a maximally compact hydrophobic core and 
hence do not crystallize in a separate transition. 

For a deeper understanding of the aggregation transi- 
tion, we now analyze entropic effects accompanying this 
transition in the microcanonical ensemble for an isolated 
system of multiple semiflexible chains with k = 0.25. The 
microcanonical entropies per monomer, s = S/N tot , are 
shown in Fig. E{a) as functions of the total energy per 
monomer, e = E/N to t- The reciprocal caloric tempera- 
ture, T- X {E) = dS/dE, is plotted in Fig. Eft). The plots 
reveal an exciting phenomenon: Increasing the energy en- 
tails a reduction of temperature in the transition region, 
known as the backbending effect [3] . This signals a phase- 
separation process [TJ|2] which is caused by surface effects 
reducing the entropy which, in an isolated system, results 
in a decrease of temperature by increasing the total system 
energy [3T11IT2"rfl"5] . This phenomenon is called "backbend- 
ing effect" , because the caloric temperature curve changes 
in the transition region its monotonic behavior with in- 
creasing total energy [3J. For truly isolated systems, this 
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2xAi 3 0.973 0.024 
3xAi 3 1.118 0.016 
4xAi 3 1.172 0.009 
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Table 1: Aggregation temperatures T agg , relative surface entropies per monomer As sur f, relative aggregation and fragmentation 
energies per monomer, e agg and ef rag , respectively, latent heat per monomer Aq, and phase-separation entropy per monomer 

Aq/T agg . 
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Fig. 3: (a) Microcanonical entropy per monomer s(e) and the 
Gibbs constructions h s (E) (dashed lines), (b) reciprocal caloric 
temperatures and Maxwell constructions, and (c) relative sur- 
face entropies per monomer As sur f. 



is a physical effect which has been verified, e.g., in atomic 
clustering experiments |16j . 

In our previous studies of heteropolymer aggregates [TJ 
|2], we took technical advantage of this microcanonical 
view to identify the peptide aggregation as so-called "cou- 
pled binding-folding transition", i.e., the individual het- 
eropolymer chains refold during the binding process and 
the finally formed aggregate possesses a single compact 
hydrophobic domain. 

By closer inspection of Fig. ^[b), we find for the ho- 
moploymer system a hiercharchical substructure caused 
by these surface effects. The number of oscillations of the 
curves, increasing with system size, reveals that the aggre- 
gation transition is actually a composition of different sub- 
processes, each of which is an individual phase-separation 
process. The amplitude of these oscillations decreases with 
system size showing that these subprocesses comprise a 
smaller surface-entropic barrier [see Fig. [3jc)] . 

The 2 x A13 system exhibits a single backbending ef- 



fect as only two chains aggregate. For the three-chain 
system 3 x A13, a different scenario is apparent. In the 
higher-energy regime, first two chains stick together, and 
the formation of the three-chain globule is a separate pro- 
cess at lower energies. This hierarchical procedure con- 
tinues for larger systems as, e.g., for 4 x A13. However, 
the impact of the individual backbending effects is getting 
weaker and, in the thermodynamic limit, these effects are 
expected to disappear asymptotically, whereas the first- 
order character of this transition remains. The horizontal 
lines in Fig. EJb) are the respective Maxwell constructions 
and define the aggregation temperatures T agg . In the fol- 
lowing, we denote the leftmost (rightmost) energy where 
T (e) = T~g as e agg (ef rag ). In the entropy curves, the 
Maxwell constructions correspond to concave hulls h s (e) 
with dh s (e)/de = T~ gg (Gibbs constructions) in the tran- 
sition regime [sec Fig.^a)], where entropy is reduced due 
to surface effects (the convex region of the microcanonical 
entropy curve is sometimes called "convex intruder" [3]). 

For a quantitative analysis, we define As(e) = h s (e) — 
s(e) which is plotted in Fig. EJc). Within the transi- 
tion region (i.e., for e agg < e < ef rag ), the peak height 
As sur f = max e agg <e<e fl - ag As(e) defines the surface entropy. 
The energetic width of the phase-coexistence regime is 
the latent heat per monomer, Aq = ef rag — e agg = 
2agg[s( e &ag) ~ s ( e agg)])- Thus the entropic phase sepa- 
ration barrier Aq/T agg should survive in the thermody- 
namic limit, if the aggregation of scmificxible polymers is 
a first-ordcr-like phase-separation process with coexistence 
of aggregates and fragments. 

Values for these quantities are listed in Table [1] for 
the three polymer systems considered in our study. We 
find that with increasing system size the surface entropies 
As sur f decrease and thus the influence of surface effects 
is getting weaker. On the other hand, the latent heat per 
monomer increases, supporting the first-order character of 
the aggregation transition. 

We have shown in this Letter that the aggregation of 
interacting scmificxible polymers is a first-order phase- 
separation process. For two reasons, we focused on small 
systems with up to four chains with 13 monomers in a 
cubic box of edge length L = 60. The first reason is that 
one of the primary goals of this study was the unraveling 
of underlying mechanisms in multiple-chain aggregation 
processes. This also includes the identification of sequen- 
tially ordered, i.e., hierarchical, subphase transitions ac- 
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companying the phase separation process of aggregation. 
As our microcanonical analysis revealed, these transitions 
exhibit peculiarities compared to thermodynamic phase 
transitions, which are due to dominant entropic surface 
effects in the small systems considered here. The surface 
changes at the interface of co-existing aggregates and frag- 
ments let the microcanonical temperature decrease while 
energy increases, whereas far away from the transition re- 
gion the more intuitive behavior is observed: temperature 
increases with energy. Because of the monotonic change of 
the microcanonical temperature in the transition region, 
this phenomenon is called "backbending" . A second rea- 
son for restricting ourselves to small systems is the demand 
for a particularly high accuracy needed for the precise es- 
timation of the density of states which is the basis for the 
microcanonical analysis. The observed surface effects arc 
expected to vanish in the thermodynamic limit [2J. How- 
ever, for molecular small-scale applications based on the 
interplay of only a few molecules and the understanding of 
biomolccular aggregation processes, the thermodynamics 
of finite-size effects is relevant. 

Our precise microcanonical analysis revealed that for 
homopolymers the aggregation transition is accompanied 
by hierarchical backbending effects. We found that after 
aggregation, complexes of rather floppy semiflexible poly- 
mers behave like globules in the liquid regime and freeze 
at low temperatures in a separate process. In contrast, for 
rather stiff homopolymers, no separate freezing transition 
is observed, in coincidence with former results |10| . The 
overall transition behavior of semiflexible homopolymers 
is also different compared to the previously studied het- 
eropolymer model [215] which only differs in the sequential 
disorder of hydrophobic and polar monomers, whereas in 
the homopolymer case all monomers are treated as hy- 
drophobic. However, we found similarities in the transi- 
tion towards the formation of a single hydrophobic do- 
main, provided the stiffness of the homopolymer chain is 
sufficiently weak. 

The understanding of aggregation and crystallization 
processes of polymers is a necessary prerequisite for the 
design of technological applications in material science as, 
e.g., ordered nanoscopic structures like fibers, nanopores, 
and channels or amorphous polymeric cells with particu- 
lar electronic properties. Aggregation processes arc also 
essential in biological systems, where enzymatic and mo- 
toric action is mediated by molecular binding processes. 
Furthermore, molecular cluster formation can also cause 
disastrous diseases like Alzheimer's, which is yet another 
reason, why generic features of polymer-polymer interac- 
tions are worth being studied. 
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